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The mean field fluctuations of large atomic ensembles can behave like bosonic modes, i.e. they induce a state 
on an appropriate system of bosonic modes. The most prominent example is that, if the atomic ensemble is in a 
homogenous product state, then the mean-field fluctuations are inducing a Gaussian state on a system of bosonic 
modes. In the present paper we show that for atomic ensemble states with exponentially decaying correlations 
(e.g. with respect to the distance of atoms) the mean-field fluctuations are inducing (possibly non-Gaussian) 
states on the on a system of bosonic modes. This result is true for a general lattice of atomic systems that is 
equipped with a reasonable distance function. 
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I. Introduction 

In many body quantum systems, the mean-field theory is focusing on the average behavior of single 
constituents (see e.g. |4 , 5, 10|). Large atomic ensembles are systems of this kind. Typical global states of 
these systems have the property that the state restricted to an single atom is independent of the individual 
atom. If we average the state restricted to a single atom over all atoms of the ensemble, we obtain the 
same state as restricted to each individual atom. For asymptotically large systems, this mean-field average 
corresponds to an effective classical systems. Since one is only concerned with expectation values of 
observable of single atoms, correlations between different atoms are irrelevant for the mean-field limit. 

A kind of "first-order correction" to the classical mean-field limit are "mean-field fluctuations". Observ- 
ables which for testing mean-field fluctuations are build as follows: One looks at the deviation of a single 
atom observable from its mean-field expectation value. Whereas the mean-field expectation value is the 
same for each single atom, the expectation value of the deviation may depend on the individual atom. A 
"mean-field fluctuation observable" (fluctuation operator) is an appropriate average of the individual mean- 
field deviations over all atoms. 

It is well known, that if a large atomic ensemble is prepared in a homogenous product state, i.e. each 
single atom is individually prepared in the same state, the mean-field fluctuations effectively behave like a 
systems of non-interacting bosonic modes. In other words, a homogenous product state of a large atomic 
ensemble, induce (via mean-field fluctuations) a "Gaussian state" on a system of bosonic modes. This 
statement has to be interpreted in the limit of infinitely large systems. This is also related to the well known 
"Holstein-Primakoff transformation" O which relates large spin systems to bosonic systems. 

The "bosonic nature" of mean-field fluctuations for a large atomic ensemble can also be interpreted as 
"simulating" bosonic systems by large atomic ensembles. 

It can be observed in experiments that mean-field fluctuations can have a "bosonic behavior" by building 
interfaces between atomic ensembles and light f2\. Here, a laser is appropriately interacting with a gas of 
atoms confined to a glass box at room temperature. The state of the laser field can be stored into the atomic 
ensemble by using the effective degrees of freedom of the mean-field fluctuations. Conversely, one can also 
perform an inverse process, by transferring the state of the mean-field fluctuations of an atomic ensemble 
to bosonic modes of a light field. 

One can also imagine to use a similar technique in order to store the state of a laser field into an ensemble 
of atoms that are trapped by the periodic potential of an optical lattice. The interesting point is here, that 
the implementation of quantum cellular automata for optical lattice systems is a natural task. In an optical 
lattice, each atoms occupies a lattice site. A quantum cellular automaton is a global (mostly reversible) 
quantum operation whose local action on a single site subsystem only affects the neighboring sites |t9J . 
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A very interesting issue combines the simulation of bosonic modes by atomic ensembles, on one hand, 
with the implementation of quantum cellular automata which acts on atomic ensembles, on the other hand: 
First, store the state of a laser field into an ensemble of atoms (arrange within a optical lattice), then imple- 
ment an quantum cellular automaton acting on the atomic ensemble, and finally "release the light" from the 
atomic ensemble. By the overall process, we obtain an incoming light field and a outgoing scattered light 
field. The question is now: 

What kind of effective operation is describing this "scattering process"? 

We expect to obtain operations that go beyond the "Gaussian world" which can be used as "non-Gaussian 
addons". Since Gaussian operations are limited in their ability to perform quantum information tasks, this 
may open a door to perform new tasks. 

Whether the mean-field fluctuations of a large atomic ensemble behave like bosonic modes depends on 
the state of the atomic ensemble. If we want to find an answer to the question given above, we will answer 
the following question first: 

For which states of large atomic ensembles do the mean-field fluctuations behave like bosonic 
modes and is there a set of states with bosonic mean field fluctuations which is invariant under 
application of quantum cellular automata? 

We shall see that Theorem [7] provides an answer to this question. For states of large atomic ensembles 
whose correlations are exponentially decaying with the distance of the single atoms (exponential clustering), 
the mean-field fluctuations behave like bosonic modes. In particular, states with exponential clustering are 
invariant under quantum cellular automata since the action of a quantum cellular automaton on a single site 
system only affects a finite set of neighbors. 

As a consequence, the following process is possible: A laser field is interacting with a large atomic 
ensemble such that the Gaussian state of the laser field is encoded into a homogenous product state of 
the atoms. A quantum cellular automaton acting on the atomic ensemble is implemented. The resulting 
state of the atomic ensemble possesses again bosonic mean-field fluctuations. Finally, we can use again 
the interaction between the laser field and the atomic ensemble to transfer the state of the atomic ensemble 
(almost faithfully) to the bosonic modes of the laser. 

The total process induces an operation on bosonic modes which maps an initial Gaussian state to some 
bosonic state, which can be non-Gaussian. With help of Theorem|5] the correlation functions of the resulting 
state can be written as a perturbation of a Gaussian state. This may be helpful in oder to decide which states 
of atomic ensembles correspond to Gaussian states. 



Outline of the paper 



In Section|ll]we provide the appropriate mathematical tools for describing mean-field fluctuation. Fixing 
a given atomic ensemble, it depends on the state which kind of system of bosonic modes (if there is one) 
is corresponding to the mean-field fluctuations. This requires to compare different bosonic systems even if 
they differ by its canonical commutator relations. The tensor algebra provides a universal description that 
covers all different bosonic systems at once. Which bosonic system is realized is part of the state of this 
"universal continuous variable system". 

How to describe and analyze mean-field fluctuation by using the framework of universal continuous vari- 
able systems (tensor algebras) is discussed in Section III Here, we also present the main results (Theorem|5] 
Theorem [7]). Technical supplements in order to give self-contained proofs are postponed to the sections in 
the appendix. 
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n. Universal description of continuous variable systems 



Our goal is to relate systems of atomic ensembles with bosonic systems. In this section, we explain this 
relation in mathematical detail and generality, by using the algebraic approach to quantum mechanics. Here 
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systems are given in terms of their observable algebras, which, in our case, are C* -algebras or more general 
*-algebras when unbounded operators are included. 

A. The tensor algebra a universal playground 

We now introduce a formalism for describing general continuous variable systems in a uniform manner, 
which is not so frequently used, but which has the advantage that the Holstein-Primakov transformation can 
be implemented easily and naturally. 

Let V he a complex vector space with a complex conjugation J. The tensor algebra over {V, J) is the 
unital associative *-algebra that is given by the complex vector space 

^{V, J) = V'®" . (1) 
The product, which is just given by the tensor product, is determined by 

(ui (g) W2 (X" • • • (8i Vn){'Wi (g) W2 ® • • • ® Wm) = t^l (g ® ' ' ' ® '^n ® W'l ® ® ' ' ' ® (2) 

and the adjoint is determined by 

(Wl (g) W2 (X" • • • (Xl Vn)* = JVn (g JUji-l ® • • • (g Jvi , (3) 

where vi, - ■ ■ ,Vn,wi, - ■ ■ , Wm € V. Note that = C corresponds to the multiples of the unit operator 
1. Obviously, there is a linear embedding $ of y into the tensor algebra J) suchthat$(w)* = ^{Jv). 
In the following, we call the operators <I>(u) "generalized field operators". 

The tensor algebra ,5^(y, J) represents the observable algebra for a wider class of continuous variable 
systems. The detailed type of the system, e.g. fermionic or bosonic, is encoded in the states under consider- 
ation. A state oj is described by normalized positive linear functional on the tensor algebra, i.e. ui{A*A) > 
and w(l) = 1. Each state is determined by the n-point correlation functions 

i^nin,- ■ ■ ,Vn) = Oj{vi (g) ■ ■ ■ (g}Vn) (4) 

where a;„ is a n-multi-linear functional on V. 

How is all this related to the ordinary Hilbert space formalism of quantum mechanics? Well, with help 
of the so caUed GNS representation we obtain a Hilbert space Hu, a vector fi^^ and a * -representation 
by linear (but not necessarily bounded) operators on with u){A) = {Q^,TTi^{A)ni^). This is just a 
consequence of the positivity of the functional oj. 

B. Realizing bosonic systems 

As a first example, let us have a look at the bosonic systems. For this purpose we construct a"quasi-free 
bosonic state" from a bilinear form, called the covariance W, on V. In order to obtain a positive functional, 
the positivity condition W{Jv, v) > has to be fulfilled. The corresponding quasi-free state is determined 
according to the following conditions: For n > we put 

wvy(ui(g)---(g)t;„) := ^ Yi (5) 
Pen2(n) {i,j)eP 

and a;^(l) = 1, where 112 (») is the set of ordered partitions of {1, • • • ,n} into two-elementary subsets. 
Note that the sum is empty for odd n. 

We associate to W the hermitian form 7 which is given by 7 (wi,t;2) = W{Jvi,V2) — W{v2,Jvi). Using 
the Araki's self-dual formalism, the self-dual CCR algebra CCR(y, J, 7) is the *-algebra that is constructed 
as follows: Let be the two sided ideal in ^{V, J) that is generated by the operators ^{v)*^{v') — 
^{v')^{v)* — 7(w, v')l. Then the corresponding self-dual CCR algebra is given by the quotient *-algebra 

CCR(y,J,7):=^(l^,J)/A- (6) 

As we will briefly sketch below, the state uw annihilates the ideal which implies that cow induces a 

unique state on CCIi{V, J, 7). Two quasi-free states ojw and ljw' on the tensor algebra belong to the same 
Bosonic system if W{Jv, v') — W{v', Jv) — W'(Jv, v') — W'{v' , Jv) = ^{v, v'). In this case, both states 
annihilates the ideal and can be lifted to the same CCR algebra. 
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Remark 1 The quasi-free states low have the special property to be "even", i.e. the expectation value of 
a single generalized field operator is vanishing ujwi^'^iv)) — 0. To obtain all quasi-free states, we take 
advantage of the following fact; Let be the real vector space of real continuous linear functionals on V . 
Note that a functional u e is real if it fulfills the condition u( Jv) — u{v) for all v ^ V. If we regard 
with its addition as an Abelian group, then is acting by *-automorphisms on the tensor algebra =^(V, J). 
For each u e V^, we define the *-automorphism q:„ G Aut{£^{V, J)) according to 

a^^iv) := ^{v) + u{v)l . (7) 

By construction, the group law is fulfilled, i.e. Q!„ia„2 = a„i+„2 is valid for all ■ui,U2 € Vg. To obtain 
a quasi-free state with a non-vanishing one-point function, we just "shift" an even quasi-free state uiw by 
an appropriate automorphism yielding the quasi-free state U!w,u = i^w ° O-u which has the one-point 
function lo\y,u{^{v)) = ujw{^{v)) + u{v) = u{v). 

C. Ideals to specify more detailed systems 

The discussion of the previous subsection shows that the tensor algebra can be indeed used to describe 
various systems by one unified object. To specify a more particular sub-class of systems additional algebraic 
relations has to be respected. This corresponds to a proper two-sided ideal ^ C ^{V,J). By inclusion, 
the set of two-sided ideals is partially ordered. As larger the ideal, as more specific is the systems class under 
consideration. For instance, if the hermitian form 7 is non-degenerate, then the ideal which describes 
the corresponding CCR-relations is maximal: This can be interpreted as the most specific description of 
a system, here for a set of bosonic modes. Each state w is accompanied with a natural system that is 
given by the quotient algebra :— ^{V, J)/ where is the two-sided ideal := {A|Vi3, C : 
u{B*AC) = 0}. Note that, by construction, does not contain the identity operator and is therefore a 
proper ideal. 

D. Comparison of states 

But what does it mean, that two states on the tensor algebra are close to each other? To give a precise 
answer to this question, we need to compare states quantitatively. For this purpose, we assume that V is 
a Banach space. The dual space of the tensor algebra ^{V, J) is denoted by ^{V, J)*. It consists of all 
linear functionals F : ^{V, J) ^ C for which for all n € N the semi-norms 

VniF) := sup |F(wi(g)---(g)t;„)| <oo (8) 

are bounded, where Vi = {v G V\ \\v\\ — 1} is the unit sphere. Now, ^{V, J)* is closed in the following 
topologies: 

• The strong topology is the locaUy convex topology that is induced by the family of semi-norms Vn, 

n eN. 

• The weak topology is the locally convex topology that is induced by the family of semi-norms 

'^{vi,- ,Vn)iP) ■= \F{vil»---l»Vn)\ with {Vi,-- - ,Vn) € Ufc 

From an experimental perspective, the strong topology is related to the comparison of two states w and w'. 
Suppose we estimate for a finite family of vectors ui , • • • , w„ £V the correlation functions io{vi^- ■ -^Vn) 
and (g) • • • (g) t)„). Then the modulus of the difference of the correlation functions can be used as a 
measure how "close" uj and uj' are to each other. 

To give an example, we consider the correlation functions of two quasi- free states low and lvw', where 
the covariances W, W have a norm difference that is given by — W\\ = sup^^ v2eVi \ W{vi,V2) — 
W'{v^,V2)\. 

Proposition 2 Let low and low' be quasi-free states with covariances W and W respectively, then for 
each n £N the semi-norm difference of the quasi-free states satisfies the bound 

n/2 

Un{u:w ^ ^w) <\\W -W'\\\n2{n)\Y, Wf-^WWr'^-K (9) 

fe=i 
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A direct consequence of the proposition (wliich we prove in the appendix) is that, if W ^ W are 
converging in norm, then low uiw' converges in the strong topology. In other words the mapping 
W ojw is continuous in the respective topologies. 

III. Mean-field fluctuations 

Many systems under consideration possessing a large number of independent degrees of freedom such 
that they can be ideahzed by infinite systems in the thermodynamic limit. Here we model this situation by 
an infinite (countable) lattice A that possesses a distance function d : A"^ ^ M+. The observable algebra of 
the global system is the so caUed quasi-local algebra that is constructed by the infinite tensor product 

2t(A) = (g) 21(A) (10) 

of single cell C*-aIgebras 21 = 2t(a;). For a given lattice point x E X, the natural embedding of the single 
site algebra which identifies 21 with 2t(x) c 21(A) is denoted by l^. We are going to use this mapping later 
on quite often. 

However, in view of non-equilibrium thermodynamics, the nature of global states of an infinite systems 
can be very different from equilibrium states and the calculation of expectation values for such states may 
be a hard computational task. The analysis of global states, one looks at the asymptotic behavior of cer- 
tain properties within the mesoscopic range. For this purpose, one restricts the global state to the local 
observable algebras that correspond to finite subsets sets X C A which is given by the finite tensor product 

2l(X) = (g) 2l(a;) . (11) 

Note that for an inclusion X C F C A, it follows immediately that 2l(X) C 2l(y). 

Taking a global state oja, we obtain for each finite subset X C A a restricted state ujx ■— i^a\%(x) 
which lives on finitely many degrees of freedom. This yields a net of sates {ojx)xca that is indexed by the 
partially ordered set of finite subsets of the lattice A. Roughly speaking, the basic idea behind the Holstein- 
Primakov transformation is to analyze the behavior of each of the states uJx concerning their "bosonic 
nature", i.e. to what extend they "simulate" continuous variable systems. We shall see, that each restriction 
LOx induces a state u)x on the tensor algebra ^(21, *), where 21 is the observable algebra of a single cell 
system. To be of "bosonic nature", the induced state uix has to fulfill "almost" the canonical commutation 
relations. This means that there is an antisymmetric hermitian form 7 such that the induced state cox is 
"almost" annihilating the ideal ^-y-. A typical behavior is uJx{A) = 0(|X|~^/^) for each operator A that 
belongs to the ideal 

A. Inducing states and y'n-fluctuations 

Let OJA be a state of the global system. Then we obtain the net of restricted states {uix)x(zA that are 
indexed by the partially ordered set of finite subsets X c A. The induction of states works by using 
"fluctuation operators" associated with the restricted state uJx and an operator a € 21: 

^ujxia) ■■= [''^'^ ~ u;x{ixa)l\ . (12) 

' ' xex 

This yields a representation $(a) 1-^ {(^) of the tensor algebra and the induced state ojx is determined 
by its n-point functions according to 

i^x{ai (g) • • • a„) := ujx{^uix{ai) ■ ■ ■ *t^x(an)) • (13) 

The main goal is now to study the asymptotic limit of large systems. For this purpose, let A be a countable 
lattice. Let {u}x)x(ia be a net of states that is indexed by finite subsets of A, where oJx is a state on 2t(X). 
The asymptotic properties in the limit X ^ A can be investigated by looking at the net of induced states 
{i^x)xcA according to the classification: 



On the bosonic behavior of mean-field fluctuations 



6 



• The state oja has -yn-fluctuations if the induced net (a)x)xcA converges w — limx^A = wa in 
the weak topology on .5^(21, *). 

• The state cja has strongly -yn-fluctuations if the induced net ((Dx)a'ca converges s — \imx-fA i^x = 
Lj\ in the strong topology on ^(21, *). 

• The state wa has weakly y^-fluctuations if for the induced net (wx)acA the each semi-norm 
n S N, is uniformly bounded: supj^^A '^ni'^x) < oo. 

Obviously, strongly yn-fluctuations implies y^-fluctuations implies weakly y^-fluctuations. We are 
now considering states that are "single site homogenous". These states are defined by the property that their 
restrictions to a single site is independent of the lattice point. 

B. Induced states for asymptotically large systems 

Asymptotically large atomic ensembles can be described by an infinite lattice system which is in some 
state uj\. Suppose we assume that the corresponding induced net of states {(I'x)xgA has weakly y/n- 
fluctuations. What conclusions can we draw from this property? What do we know about the asymptotic 
behavior of the correlation functions of the induced states lix? 

Since each semi-norm i^„(a)x) is uniformly bounded in the size of the subset X, we know that there are 
weak limit points. In order to analyze the properties of these limit points more systematically, we will give 
here an "operational" description of what limit points are. 

Within a concrete experimental realization, the atomic ensemble under consideration will be always 
finite. If the setup is scalable, then, at least in principle, the same experiment can be performed for various 
sizes of the system, i.e. the subset X can be regarded as a "classical configuration". Here one can also 
think of a situation, where atoms occupy only finitely many sites of a lattice randomly. Thus we are dealing 
with a preparation device that prepares for each finite atomic ensemble X C A a state uix with a certain 
probability /i(A'). The probability distribution : X i— ii{X) is nothing else but a classical state on the 
system of finite subsets X C A. The corresponding observable algebra consists of all bounded complex 
valued functions / : A D X i— > f{X). The expectation value of an observable / for the state /i is then 

given by a*(/) = ExcaM^)/W- 

A general classical state on the system of finite subsets AT C A is a complex valued hnear functional on 

the algebra of bounded complex valued functions f : X ^ f{^) such that the following holds: 

• Positivity: r]{f) > for each / > 0. 

• NormaUzation: 77(1) = 1. 

A preparation device that produces asymptotically large atomic ensembles has the property that, in the 
limit X -> A, the probability that only a finite number of lattice sites are occupied is vanishing. This 
corresponds to classical states 77 with the following property: 

• ?7(/) -Oiflimx^A/W = 0. 

A state 1] with this property is called a "limit point". To justify this notion, suppose that limit limx f{X) = c 
exists. In this case \imx {f {X) — c) — 0, and rj{f — cl) ~ vif) — c — Q follows, which means that there 
expectation value ri{f) = c coincides for all limit points 77. 

What can we say about the limit points of the induced net {ujx)xcA for a state lja that have weakly 
y^-fluctuations? To each operator A E >5'(2t, *) of the tensor algebra, we assign a bounded function 
which is given by AT i-> u}x{A) iflTl . Now, each limit point 77 induces a state on the tensor algebra by 
uJrjiA) = rj{X H' ujx{A)). Here we use the suggestive notation rj{X 1— > f{X)) := rj{f ) to represent an 
expectation value. 

The states a),, describe the mean-field fluctuations of asymptotically large atomic ensembles. The next 
propositions states that these mean-field fluctuations behave like bosonic modes. Consider a state wa 
that is single site homogenous with single site restriction cj = uj^ ° ix- Then there is a natural an- 
tisymmetric hermitian form 7(0,6) := a;([a*,6]) on the observable algebra 21 of the single site sys- 
tem. The ideal which represents the canonical commutation relations, is generated by the operators 
/^(a,6) :- [<i>(a),<i>(6)]-7(a*,fo)l. 
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Proposition 3 Let lu\ be a single site homogenous state having weakly ^Jxi-fluctuations. Then for each limit 
point rj, the state uj^ annihilates the ideal and can uniquely he lifted to a state on the corresponding 
CCR algebra. 

Proof: By Lemmal9]of the appendix, we conclude that limx i^x (^) — for each operator in the ideal 
This implies Wr^{A)— r]{X t-^ ujx{A)) — which implies that o)^ annihilates ■ 

Remark 4 Proposition |3]can be interpreted, at least to a certain extend, by saying that a state uj\ of a large 
atomic ensembles with weakly i/n-fluctuations possess bosonic mean-field fluctuations. This is justified 
by the fact that each limit point cj^ is a state on the CCR algebra CCR(2t, *, 7) which describes a bosonic 
system. On the other hand, the CCR algebra is an algebra of unbounded operators and it might happen 
that the GNS representation associated to a limit state a),, has "exotic" properties. Recall, that the GNS 
representation is given by a Hilbert space H an algebra homomorphism tt that assigns to each operator 
A in the CCR algebra a linear (unbounded) operator on H as well as a normalized vector il £ H such 
that uJrf{A) = {n, TT{A)n). The question that arises here is whether it is possible to build the exponential 
exp(i7r($(a))) of a field operator 7r($(a)) in the representation tt, where a = a* is selfadjoint. If we can 
do this, then we obtain a representation of the Weyl algebra by bounded operators. If it is not possible to 
build the exponential we are dealing with an "exotic" case (see e.g. |6|). In order to exclude this kind of 
pathologies, we need to consider more specific examples of states. 

C. States with exponential clustering 

A state WA on 21(A) has exponential clustering (with respect to d) if for local operators A G 2l(Ar) and 
B e 2l(r) the identity 

oja{AB) - ^a(A)^a(B) + G(x,Y) (A, B)e-'^(^'^) (14) 

is valid for a bounded bilinear function G^x.y) ■ 2l(A') x 2t(F) C such that \G(x.Y)iA, B)\ < 
Go||v4||||_B|| for all A,B, for all finite regions X,Y. Here Gq is a constant that is independent of the 
localization regions. The bilinear forms Gx,y express locally the deviations from the state to a product 
state, being scaled with the exponential of the distance. Therefore G(x,y) indicates the presence of corre- 
lations that are exponentially decreasing with the distance. To give a name, we call the family of bilinear 
maps G — (G(x.y))x.YcA the correlators. Note that, equivalently, exponential clustering is given by the 
condition 

\iOA{AB) - u;a{A)u;a{B)\ < Go e-'^^^^^) (15) 

for all A € 2l(A'),B e 2t(y). Here d{X,Y) = mmrc^x,yGY d{x,y) is the distance between the finite 
subsets X,Y C A. We always require here, that the distance d is regular, i.e. the maximal number N{r) of 
lattice sites within a ball of radius r is bounded by a polynomial. 

The exponential clustering property can be used to derive a useful cluster expansion in terms of expecta- 
tion values of the single site restriction cj and the correlators G. In order to write down this expansion, we 
introduce the following objects: 

• For each finite subset y C A we introduce the "the spread" A(y) := max^gy d{y, Y \y) which 
measures the maximal distance of a point to its relative complement in Y. 

• An fc-elementary subset {yi, • • • , yk} C X is called "spread optimally enumerated" if the enumera- 
tion fulfills the condition d{yi, {yi+i, ■ ■ ■ , yk}) — ^{yi, ■ ■ ■ , yk) for alH = 1, • • • ,k — 1. Note that 
each subset can be spread optimally enumerated. 

• Given a tuple x E X", we choose a spread optimal enumeration of the range Ran(a;) = 
{yi,--- ,2/|Ran(x)|} and we consider the correlators G^ := G(^y^^{y^^^^... ^y^.^^^^^^^}) which test 
the correlations for splitting the site yk from the remaining points {y^+i, • • • ,y|Ran(x)|}' where 
fc = I,-- • , |Ran(a;)| - 1. 
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• For a family of operators oi, • • • , a„ G 21 and a tuple x G X" whose range {yi, • • • , y|Raii(x)|} 
is spread optimally enumerated, we introduce the single site "cluster operators" G 21 which are 
given by the ordered product :— Yijex-'^iyk) where the ordering is according to the value of 
the index in a;"i(?/fc) = {j 1, • • • ,n\xj ^yk}- 

The following theorem, whose proof is given in the appendix, states that correlation functions of the 
induced states Cjx admit a cluster expansion in terms of the single site restriction uj, the correlators G and 
cluster operators a^: 

Theorem 5 (Cluster expansion) Let wa be a single site homogenous state with single site restriction uj 
and exponential clustering with respect to d. For each fli, • • • , a„ G ker(a;) and for each finite subset 
X C A the n-point correlation function of the induced state ux can be written as 

ujx{ai® ■ ■ ■ ® an) = io'^^iai (g) ■ ■ ■ (g) a^) + Fx{ai (g) • • • (8 a„) . (16) 

where the correlation function of the induces homogenous product state ui'^^ and the functional Fx are 
given by 



UJ 



|Ran(a:)|-l 

Fx(ai0---®a„) = iXj-t ^ ^ uj{al)---uj{al,,) ^^^^ 

icex" k=i 

where for each x G X" the range Ran(a;) is spread optimally enumerated. 

As the cluster expansion is stated above, it holds for all correlation functions for which the operators 
ai, • • • ,a„ are chosen in the kernel ker(u;) of the single site restriction. If this is not the case, we can 
express the correlation function in terms of ai = a[ + w(ai)l where a'^ G ker(w). The tensor product 
ai (8) • • • (8 a„ can be expanded in terms of the operators a'^ G ker(a;) according to 

ai (8 •••«)«„= ^ 0a- W u}{aj) (18) 

,/C{l,--- ,n} ieJ ,n}\J 

where the sum runs over all ordered subsets. To get the general cluster expansion for the full tensor algebra, 
one only has to apply Theorem|5]to (18 1 for each summand. 



It is known, that the induced net Juj®'^)xcA converges weakly to a quasi-free state. We show here a 
slightly stronger result; 

Proposition 6 A homogenous product state uj^-^ has strongly ^/n-fluctuations. In particular, the induced 
net {Cj'^'^)x(ih converges strongly to the quasi-free state uj^f whose covariance is given by the truncated 
two-point function W{a, b) = uj{ab) — a;(a)w(6). 

Homogenous product states are the simplest among states that have exponential clustering. For the 
general case, the following is true: 

Theorem 7 Each single site homogenous state with exponential clustering has weakly ^/ri-fiuctuations. 

The proof of the theorem is quite technical and therefore postponed to the appendix. However, it takes 
advantage of the cluster expansion of Fx into single site expectation values and correlators. The basic idea 
to get a uniform bound for the semi norms {Fx ) is to count the number of terms that are contributing to 
the cluster expansion. In total, we sum over all tuples in A"" which gives jATj" terms. Since we normal- 
ize by multiplying jATj^"/'^, a naive counting would give the non- uniform bound z^„(Fjf) < 0(|X|-"/2). 
By a more careful analysis, it turns out that effectively only |A'|"/^ terms are contributing. By choos- 
ing ai,-- • , a„ G ker(aj), the single site expectation value of a cluster operator w(a^) is vanishing if 
x~^{yk) = {j} contains only a single element. Note that in this case we just have w(a^) = io{aj) = 0. 



This reduces directly the number of terms which in the cluster expansion ( 17 1. A large number contributions 
are also coming from tuples x with range {yi, ■■■ , ?;|Ran(i;)|} for which the spreads A{yk, • • • , y|Ran(x)|) 
are large. These contributions are also of order jATj"/^, since they are suppressed the exponential damping 
exp(-A(yfe,--- ,y|Ran(x)|))- 
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Remark 8 We can derive from the cluster expansion that for asymptotically large atomic ensembles there 
is a quasi-free part from the product state contribution and a perturbation which comes from the correlators. 
Namely, for each weak limit point ry the state can be written as 

LO,^ = OJql + F,, (19) 

with F,i{A) = ri{X i-> Fx {A}). The functional F,, is a perturbation of the quasi-free limit state Wqt which 
may depend on the limit functional ij. Note that Theorem |7] guarantees the existence of weak limit points 
F^, since sup^cA ^'n(-Fx) < oo. 

IV. Conclusion 

We have shown that states of large atomic ensembles whose correlations are exponentially decaying with 
the distance between atoms (exponential clustering) possess bosonic mean-field fluctuations. In addition to 
that, these states are invariant under applications of quantum cellular. 

This enables the implementation of the following type of process: The bosonic modes of a light field are 
coupled to a large atomic ensemble such that the Gaussian state of the laser field is transferred almost per- 
fectly (where the precision is here of order 0(^/number of single atom systems)) to a homogenous product 
state of the atoms. A quantum cellular automaton acting on the atomic ensemble is implemented. The re- 
sulting state of the atomic ensemble possesses again bosonic mean-field fluctuations and the resulting state 
of the atomic ensemble can be transferred back almost perfectly to the bosonic modes of the light field. 

The total process induces an operation on bosonic modes which maps an initial Gaussian state to some 
bosonic state, which can be non-Gaussian. With help of Theorem|5] the correlation functions of the resulting 
state can be written as a perturbation of a Gaussian state. This may be helpful in oder to decide which states 
of atomic ensembles correspond to Gaussian states. 

It is still an open problem to decide in general from the state of the atomic ensemble whether the resulting 
induced state is Gaussian or not. Concerning states with exponential clustering, the cluster expansion 
(Theorem |5]l appears to be a reasonable technique in order to address this problem. Here the correlation 
functions of the fluctuation operators can be expanded into the correlation functions of the homogenous 
product state uj'^'^ (here ut is the restriction of the global state to a single atom) and some correction Fx- 
For large atomic ensembles, the correlation functions of the homogenous product state uj'^^ correspond to 
a Gaussian state, whereas Fx can be regarded as a "perturbation". 

Furthermore, it would be desirable to construct new examples of atomic ensemble states (in particular 
beyond homogenous product states) whose induced net has strongly y^-fluctuations or y^-fluctuations. 

In order to archive more concrete results in this direction, one has to consider here more concrete exam- 
ples. One suggestion is to consider ensembles of two-level atoms arranged in a one-dimensional lattice. A 
natural class of states for which mean-field fluctuations can be investigated are stabilizer states which are 
invariant under the action of so called Clifford quantum cellular automata (see |[il|7][8l and references given 
therein). 
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[11] This function is indeed bounded which can be verified as follows: The operator A can be written as a fi- 
nite direct sum 0^^,) Ak with Ak G 21®''. Since lja has weakly ^-fluctuations we obtain that |a;x(j4.)| < 
Y,k=o \'^x{A„)\ < X^^^o Cn||A„|| < 00 with C„ = sup;^^-^ fn(wx). 
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A. On quasi-free states (proof of Proposition |2| 

Proof of Proposition^ Let W, W be two bounded covariances (positive bounded bilinear form) on V. 
Recall that a bilinear form F on y is bounded if ||i^|| :— sup,„^ uaeVi W{'^It''^2)\ < oo. Moreover, recall 
that a valid covariance W has to be positive in the sense that W{Jv, v) > for all v E V. We fix vectors 
vi, - ■ ■ ,Vn and use the expression (jSj for the correlation functions of quasi-free states to calculate the differ- 
ence of the n-point functions for the quasi-free states ujw and tuw' ■ For each ordered partition P e 112 (n) 
we choose an enumeration P = {{ii, ji}, ■ ' ' i {*n/2) Jn/2}) and we introduce for Z = 1, • • • ,n/2 the quan- 
tities Wpj :— W{vi^,Vj^) and W'pi :— W'{vi^,Vjj). This yields for the difference of the coiTesponding 
correlation functions: 



(Al) 



UJw{vi ® • • • ® Vn) — UJw'{vi ® • ■ ■ <8) Vn) 

n/2 

Pen2(n) /=i 

Here we have used that the difference of products can be written as a sum in the following way: 

n/2 

t^p.i • • • W^p,„/2 - • • • W^;,„/2 = E ^^.1 • • ■ (W^P.; - M^^.i)M^p.i+i ■ ■ • W^p,„/2 ■ (A2) 

By using the fact that W and W are bounded bilinear forms, the modulus of Wp^i (and similarly for W 
and — W) can be bounded by | Pl/p ; | < || H'l'jJI || VK|| and we obtain the desired bound 

< \\v,\\ ■ ■ ■ \\v4 \\w- w'w \u,{n)\J2 Www'-'ww'r^'-' . 

1=1 



B. Proving the bosonic behavior of mean-fleld fluctuations for large atomic ensembles 

Roughly, the statement of Proposition |3] is that mean-field fluctuations for states having weakly y/n- 
fluctuations behave like a bosonic system for large atomic ensembles. The following lemma provides the 
bounds which are used to prove this statement. 

Lemma 9 Let uj\ be a single site homogenous state on the quasi-local 21(A) such that the induced net 
{'^x)x<iA has weakly ^Jii-fiuctuations. Then for each family of operators ai, • • • , a„ the bound 



($(ai) • • • $(a,_i)/^(a„ a,+i)$(a,+2) • • • Han))\ < 2\X\-^/^ Cn-i [] (Bl) 



holds with Cn ■= SUPjs^^^ l/n^LOx)- 

Proof: The fluctuation operators define a representation 7rtj^($(a)) := <i>„^(a) of the tensor algebra 
■9'(2t, *) by operators in 2t(X). The induced state ujx = ° t^ujx is just the pullback of the state 
ujx by the representation tt^^ . According to the definition of fluctuation operators, the identity 

^^,(/^(a,6)) = [$^^(a),<i>„,(6)]-7(a*,fo)l = \X\-'/^^^x{[a,b]) (B2) 
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is vaUd for all single site operators a, b. Inserting this identity within the correlation function for oi , • • • , a„ 
implies 

|wx(*(ai) • • • $(ai_i)/^(ai, aj+i)$(a,+2) • • • $(an))| 

= \X\-'^^\ujxiHai) ■ ■ ■ $(a,_i)$([a„ a,+i])$(a,+2) • • • $(a„))| 

<2\x\-'/^c„_^]lM\ . 

i=l 

which proves the proposition. Recall that we have used that all the semi-norms ^'„ are uniformly bounded. 



C. Cluster expansion for correlation functions (proof of Theorem|5j 

This subsection provides the proof of the expansion (Theorem |5]) of correlations functions for states 
cja with exponential clustering. We also assume here that wa is single site homogenous with single site 
restriction w. 

We first derive here an simpler expansion for correlation functions of the form LOf^{Ai • • • Aj.) — 
ujy{Ai ■ ■ ■ Ak), where Y = {yi,y2,'-- ,yk} C A is spread optimally enumerated and the operator 
Ai G is localized at site y^. The idea of the expansion is to expresses the expectation value 

ujy{Ai ■ ■ ■ Ah) in terms of single site expectation values i^(^j) as well as terms Gj{Aj, Aj^i ■ ■ ■ Ak) :— 
^({wl-iw+i.- - -Vk}) (AjjAj^i ■ ■ ■ Afc) that are given by the bilinear forms G'(x.y)- Recall, that for two 
operators A E 21(X) and B e 21(F) the exponential clustering can be expressed by the identity 

ujA{AB)=iJx{A)ojY{B) + G(^x.YM,B)e-''^'''''^ . (CI) 

We use the following lemma in order to prove Proposition|5] 

Lemma 10 Let Y — {j/i, • • • , yk} be spread optimally enumerated and let Ai G 2l(yi), i = 1^ ■ ■ ■ ^k, be 
single site operators. Then the expectation value loy [Ai ■ ■ ■ A^) can be expressed as 

ujY{Ai---Ak) = uj{Ai) ■ ■ ■ uj{Au) 

fC2) 

+ Y^u:{A,)---^{Ai_,)Gi{Ai,Ai+,---Ak)e-^^y--^y-^ 
1=1 

where the bilinear form Gi is defined as given above. 

Proof: Let Y = {yi, ■ ■ ■ , j/^} be spread optimally enumerated. Then {y2; ■ ■ ■ , j/fc} is also spread optimally 
enumerated. Suppose now the expansion is valid for spread optimally enumerated sets with k — 1 elements. 



UJY\yi (^2 • • • ^fc) = '^(^2) • • • w(Afc) 

''"^ (C3) 
+ Y,^iA2)---uj{Ai_,)GiiAi,A,+,---Ak)e-^^y"-'y''^ 

1=2 

where the product ^(^2) • • ■ ^(^i-i) = 1 is declared to be empty for / = 2. Then we can use the expansion 

Wy(^lA2 • • • Afc) - Uj{A^)^Y\y, (^2 ' ' ' ^fe) + Gi(Ai, ^2 ' ' ' ^fc )e-^(«i ^^"'^ (C4) 

By inserting the expansion for i-^Y\yx (^2 • ■ ■ ^fe) gives 
^y{Ai---Au) = w(Ai)cj(A2)---a;(Afe) 

+ ^(Ai)c^(^2) • • • ^{M^i)Gi{AuAi^^ ■ ■ ■ Ak)e-^(y'- ^y^^ 
1=2 

+ GiiA,,A2---Ak)e-^'^y''-'y-^ (C5) 

= uj{Ai)oj{A2) ■ ■ ■ uj(Ak) 
fc-1 

+ uj{A,MA2) ■ ■ ■ uj{Ai^i)GiiAi,Ai+i ■ ■ ■ Ak)e'^^y'-- '^'=) 
1=1 
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Note that for a one-elementary set the statement is trivial and for a two elementary set , j/2 } the expansion 
is also valid since we have uj\{AiA2) — a;(yli)a;(^2) + ^'(yj^j,^) (Ai, yl2)e~'''^i'^^^ and for any two 
elementary set the distance coincides with the spread. ■ 



Proof of Theorem^ Let uja be a single site homogenous state with strong exponential clustering. Then for 
operators ai, 02, • • • , a„ S ker(w) the correlation function for the induced state ujx can be written as 

cix(ai ® ■ ■ ■ ® a„) = ^ UJX (<fl2 ■ • • afRa„(x)|) • (C6) 

where af. — Yijex-'^iv ) cluster operators. If we enumerate for each x E X" the range Ran(x) 



spread optimally, we can apply Lemma 10 to the expectation values 



UJX (0.10.2 ■ • ■afRan(x)l) = '^(^1) ' ' ' '^(afRan(x)| ) 



|R,an(2;)|-l 

+ Yl ^K)---uj{al^i) ^C7) 

k=l 

X Glial, al^, ■ ■ ■ afR,„(,)|)e-^('"=-'"=+-- '^i-"<^)i) 

where Gl = G({j^^}.{y^^j_... The statement of Proposition jSjfoflows directly by summing over 

the elements in X" and normalizing by ■ 



D. Remarks on the strong topology of the tensor algebra 

Let cj be a state of a C*-algebra 2t. Besides the strong topology, we introduce here the "w-strong topol- 
ogy" on the space of continuous linear functionals 5^(2t, *)*. It is defined to be induced by the family of 
semi-norms v!^, n E N, where i'!^ assign to each functional F £ ^(21, *)* the value 

iy!:l{F):= sup ||ai|r'---||a„irV(ai® . (Dl) 

ai ,aTifEkcr(a;) 

The reason for introducing the semi-norms i/^^ is that for the functionals we are dealing with (correlation 
functions of fluctuation operators) the semi-norms are easier to estimate. In view of this, the following 
lemma is helpful: 

Lemma 11 For a given state uj on 21, the uj-strong topology is equivalent to the strong topology on 
,^(21, *)*. In particular the bounds 

<<^n<Y. f 2^= (D2) 



k=0 

are valid for all n g N. 



Proof: Each operator a G 21 can be written as a' + cl with a' £ ker(a;) and c = uj{a). This yields 

ai (g) ■ ■ ■ (E) a,i = {a[ + cil) (g) ■ ■ ■ (E) {a^ + Cnl) ^ ^ H • 

IC{l,2,---n} jei" iei 

Applying the linear functional on both sides and taking the modulus, we obtain 

iF(ai®...®a„)i < n ii^^'unii^^ii^/K^) 

7C{l,2,---n} je/'- iei 

<||ai||...||a„|| J2 2l^'</|(^)- 

/C{l,2,-Ti} 



(D3) 



(D4) 
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Here we have used that ||a^|| < 2||ai||. This implies the desired bound 

7C{l,2,-n} 

Note that the semi-norm i/^ is optimizing the modulus of the value of a functional over tuples (ai, • • • , a„) 
in the kernel ker(a;). This immediately implies that i/„ > i^!^. ■ 

E. States with exponential clustering have weakly -/n-fluctuations (proof of Theorem [tJ 
1. Homogenous product states (proof of Proposition |6| 

Proof of Proposition^ Let uj\ — u!^^ be a homogenous product state. For operators ai, • • • , a„ G ker(aj) 
the correlation function of the induced state ujx is given by 

w®^(ai(^---(^a„) = 51 ^(°?)-"^KRan(.)|) > (El) 

where is the cluster operator = Y[jex-'^{yk) ^■i- Ek{X) the set of all enumerated fc-elementary 
subsets in X and let n(fc, n) the set of all ordered partitions of {1, • • • , n} into k non-empty subsets. Then 
( |E2[ ) can be written as 

n 

uj®^{ai®---®an) = \X\-'^Yl '^(a/i)---^(«/J 

{yu- ,yk}eEk{x) (/i,- ,/fc)en(fe,n) 

k=i (/i,--- jfc)en(fc.n) 

with cluster operators ai- := Ilje/ '^^^ Here we have used the fact that for each tuple x E X" there is 
a unique enumerated fc-elementary subset {yi, • • • , i/k} (k < n) and a partition (/i, • • • , Ik) G n(fc, n) 
such that Xi — yj for i <E /j. Suppose that for a partition (/i, • • • , /fc) one of the sets Ij contains only one 
element Ij — {/}, then this partition is not contributing to the sum since uj{ai- ) = ijj{ai) = 0. Hence, only 
those partitions with \Ij \ > 2 are contributing. For the n-point correlation function (restricted to ker(a;)'^") 
we can write 

1=0 ^ 11^ (/i,-,/„/2)en2(n) 



\X,, 

k=l,k^n/2 1=0 ^ ' (/l,--- ■-ffc)en>2(fe,n) 



(E3) 



where 112 (n) is the set of all ordered partitions into two elementary subsets and n>2(fc, n) is the set of all 
ordered partitions into k subsets containing more than one element (/i, • • • , Ik), where at least one subset 



Ij contains more than two elements. Note, that the first term in (E3 1 is vanishing if n is odd. Ek{X) is the 

set of enumerated fc-elementary subsets whose cardinality is = nf=o^(l^l ~ 0- Restricted to the 

sub-algebra ^(ker(a;), *), the functional oj'^^ can be written as 



V2-1 , ^ 



(=0 



with a continuous functional Dx ^ ^(21, that is given by 

Dx(ai®...®a„):=^|Xr"/2+^-n (l-^) E uj{aj,) ■ ■ ■ uj{ajj (E5) 

k=l 1=0^ 11/ (/i,... ,/fc)gn>2(/£,n) 
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for ai, • • • , a„ € ker(cj). To bound the semi-norm v':^{Dx), we observe the bound 

n k — 1 



\Dx{,ai®---®an)\ < llaill •••||a„|| ^ |X| 



-ri/2+fe 



k=l 



1=0 



I 

\X\ 



|n>2(fc,n)|. (E6) 



If we consider a partition (/i, • • • , Ik) in n>2(fc, n) then the constraint 2fc + 1 < n has to be fulfilled. This 
implies < |x|-i/2 and we obtain: 



n k — 1 



\Dx{ai (g) • • • a„)| < llaill • • • ||a„|| \X\-^/^ II 11 ( ^ " m ) ''^'^^^^ + 1 < ") ' (E7) 



fc=i ;=o 



where n) = |n(fc, n)| is the number of all partitions of the set {1, • • • ,n} into k non-empty subsets 
(Stirling number of second kind). For a logical statement S, the (5-function is defined by 5{S) = 1, if S' is 
true, and 5{S) = i f S i s false. By Lemma 11 it is sufficient to show convergence for the semi-norms j/^. 
From the inequality (E7 1 we obtain for |X| > n the bound 

n 

<{Dx)< |Xri/2^5(fc,n)5(2fc + l <n) 



(E8) 



fc=i 



which implies that s — lim^cA Dx = 0. Since limxcA Iir^o ^ ^ jX| ) ^ proposition follows 



finally from (E4 1: 



s - lim uj^^ 

XCA 



n/2-1 



lim TT 1 



1=0 



I 

\X\ 



ri/2-1 

lim TT (l 



I 

m 



l^qf = Wqf . (E9) 



2. Estimating the number of subsets for a given spread 

To formulate our next lemma, we introduce the function A that assigns to each finite subset Y C A the 
maximal distance of a point in Y to its complement: A{Y) := maxy^y d{y, Y \ y). To prove bound on 
correlation functions, one partial task is to count for a finite subset X C A the number N{X, k, r) of all 
/c-elementary (k > 2) subsets Y in X such that A(F) < r. Whereas the spread A measures the spreading 
of one point sets, we can also look at the spreading of sets that contain more elements. We introduce 
the fc-spread of a set Y as Afe(y) = maxjgpj,(-y) d{J, Y \ J) which measures the largest distance of a 
A;-elementary subset J £ Pk{Y) to its relative complement in Y (here Pk{Y) denotes the set of all k- 
elementary subsets in Y). According to this definition, the one-spread Ai = A is just the spread. 

Lemma 12 Let Y be a subset of X with A{Y) < r and IS.2{Y) > r. Then Y = Z U {x, y} is the disjoint 
union of a set Z and a two elementary subset {a;, y} such that A(Z) < r and d{x, y) < r. 



Proof: Since the 2-spread of Y is larger than r, there exists a pair of points {x, y} such that d{{x, y},Y \ 
{x, y}) > r. Let Z = Y \ {x, y} be the relative complement of {x, y}. Then the spread of Y can be 
expressed as 

A(y) — max {min{d(2;, Z \ z), d{z, x), d(z, y)}, miii{d{x, Z), d{x, y)}, mm{d{y, Z), d{x, y)}} . 

(ElO) 

Since A(y) < r the bound mm{d{y, Z), d{x, y)} < r has to be fulfilled. This implies that d{x, y) < r 
because d{y, Z) > r holds by our assumption. Moreover, the inequality 

A(r) > max{min{d(z, Z\ z),d(z,a;),d(z,y)}} (Ell) 
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is obviously fulfilled. For each z e Z we have r > min{c?(z, Z \ z), d{z, x),d{z, y)} ~ d{z, Z\z) since 
d{z^ x) > r and y) > r is valid according to our assumption that each point in Z has a distance > r to 
X and y. This implies the inequality 

r> A(y) >maxd(z,Z\z) = A(Z) . (E12) 

zeZ 

Hence, Y can be decomposed into a disjoint union of a set Z and two points {x, y} such that A(Z) < r 
and d(a::, y) < r. ■ 

Lemma 13 Lef Y be a subset such that A(y) < r and IS.2{Y) < r. Then Y — Z U {x} is the disjoint 
union of a set Z and a single point x such that A(Z) < r. 

Proof: Let A2{Y) < r, i.e. for each pair {x,y}, the distance to the relative complement d{{x,y},Y \ 
{x, y}) < r. In order to discuss this case, we introduce for each point y E Y the number I{y, Y,r) — \{z 
^ \ y\d{y, z) < r}\ of points in the relative complement of y whose distance to y is smaller than r. Now 
we have to perform a further case distinction: 

• There exists a point x E Y such that I{x, Y, r) — 1. In this case the spread of Y can be expressed in 
terms of the set Z = Y \ x and the single point x as 

A{Y) = max {mm{d{z, Z \ z),d{z,x)},d(x, Z)} . (E13) 

Suppose now that d{x, z) < r holds. Then this is true for a unique point z — z* £ Z. Since the 
2-spread fulfills the bound A.2(Y) < r, for all y S Z\ z* the bound d{{x, z*}, y) < r follows. This 
implies that d{z* , y) < r since d{x, y) > r holds for all points y £ Z\z* . Hence d{z* , Z \ z*) < r. 
In all other cases, we have d{x, z) > r which implies r > mm{d{z, Z \ z), d{z, x)} — d{z, Z \ z). 
Putting these things together implies finally 

A(Z) — max(i(z, Z \ z) < r . (E14) 

• For all points y e F we have /(y, Y, r) > 2. In this case we can take any point x and express the 
spread of Y as 

A{Y) —nmx{mm{d{z, Z \ z),d{z,x)},d(x, Z)} . (E15) 

with Z = Y \ X. Suppose now that d{z, x) < r for some z E Z. Since I{x, Y,r) > 2 there exists 
at least one y E Z such that d{z, y) < r which implies d{z, Z \ z) < r. On the other hand, if 
d{z, x) > r, then d{z, Z\ z) < r since mm{d{z, Z \ z),d{z, x)} < r. This implies again 

A{Z) = niaxd(z, Z\z)<r . (E16) 



Lemma 14 For each finite subset X,for each r > 0, and for each n > 2 the bounds 

N{X,n,r)<q[,,]\Xr/^N{rr^^ (E17) 

are valid, where the numbers are recursively determined by q[n+2] = {n + \)qyn+i\ + with initial 
conditions gp] = 1 and gpj = 2. 

Proof: According to Lemma [T2| and Lemma [13] we have the following two cases: 

1. Each subset y C X with A{Y) < r is a disjoint union Y = Z \J {x, y} with A{Z) < r and 
d{x,y) < r. 



2. y is a disjoint union Y = ZU {x} with A{Z) < 



17 



On the bosonic behavior of mean-field fluctuations 



Putting both cases together, we obtain a recursion formula for bounding N{X, n + 2,r). For the first case, 
we can choose for each n-elementary subset Z with A(Z) < r a pair {x,y} with d{x,y) < r. This 
gives less than N{X, n, r)\X\N{T) possibilities. We have to add the second case, where for each n + 1- 
elementary subset Z with A(Z) < r we can just add a point x such that A{ZU {x}) < r. Since there must 
be a point z E Z with d{x, z) < r, we have (n + l)iV(r) possibilities to choose x which gives less than 
{n + l)N{r)N(X, n + l,r) possibilities. In total, we obtain the recursive bound 

N{X, n + 2, r) < {n + l)N{r)N{X, n + l,r) + \X\N{r)N{X, n, r) . (E18) 

To get a convenient explicit solution of this recursion, we allow to over-count a bit by using the inequality 

N{r) < A^(r)i/2|x|i/2. 

Nix, n + 2,r)<{n + l)N{ry/'^\X\^/'^ N{X, n + 1, r) + \X\N{r)NiX, n, r) . (E19) 

Now, we insert the ansatz N{X,n,r) < |X|"/^A^(r)"/^ into the recursion bound which gives the 
following consistency relation: 

g[„+2]|^|^"+'^/'iV(r)("+')/' < (n+ l)Ar(r)i/2|X|i/\[„+i]|X|("+i)/2iV(r)("+i)/2 

+ |X|7V(r)(7[„]|X|"/27V(^)«/2 (E20) 
< + +?[„]] |X|("+2)/27V(r)("+2)/2 

Thus, if the sequence (/[„] fulfills the recursion relation qin+2] = {n + l)q[n+i] + qin], then we obtain a 
consistent upper bound. For getting the correct initial conditions with respect to n, we first look at the case 
k = 2. The first point of the requested set can be chosen freely in X, which gives \X\ possibilities. The 
second point must have distance r to the one firstly chosen. Hence we obtain the bound N{X, 2, r) < 
|X|A^(r). For fc = 3 two points within the requested set must have distance r. This gives at most |X|iV(r) 
possibilities. The third point must have a distance < r to one of the chosen two points. This yields 
N{X, 3, r) < 2\X\N{ry < 2\X\^/^N{r f/^. ■ 



3. Towards bounding the semi-norms Vn 

Recall that A is a countably infinite lattice with a regular distance d. Here regular means that the maximal 
number of lattice sites within a ball of radius r is bounded by a polynomial P{r), i.e. N{r) = sup^.^^ |{y e 
A\d{x, y) < r}\ < P[r) for aU r e M+. 

If we require that the distance d is regular, than N{r) is bounded by a polynomial which implies that 
Nk ~ YlT=i N{r)^e~^ is finite and monotonous increasing in k. For each finite subset X, the semi-norms 
Vn of the induced state Cjx = cj®"^ + Fx can be bounded bybounding the corresponding semi-norms of 
Qj^x ^jjj separately. We know already from Proposition |6] that the semi-norms of Cj®-^ are uniformly 
bounded in X. Thus it is sufficient to bound the semi -norms of Fx- For this purpose, we introduce the 
quantities 

|R,an(i;)|-l k-1 
ajgjf" k=l 1=1 

Lemma 15 For each finite subset X d A, Bn{X) fulfill the bound 

Bn{X) < B„ |Xr/2 (E22) 
for each finite subset X d A, where are finite positive numbers that are given by 

71 k — 1 

Br^^nlY, S{k, n) ^ q[k-i+i]Nk-i+i . (E23) 

k=l 1=1 

Here S{k, n) are the Stirling numbers of the second kind and the numbers are recursively determined 
by q[n+2] — {n -\- l)(7[„+i] + q^^] with initial conditions gp] = 1 "^n^^ 9(3] = 2. 
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Proof: To bound Bn{X), we take advantage of the following fact: Let P{n,k) be the set of partitions 
of n into k non-vanishing summands. For each tuple x G X" there exists a fc-elementary subset Y — 
{Ui, • • • , Uk} C X with k < n and a partition (ni, • • • , Uk) E P{k, n) as well as a permutation a E Sn 
such that X = <T(?/f , • • • , 2/^"*° ) where y^' is the /-tuple with constant entry y. Note that the permutation 
(7 is not unique, since permuting one of the sub-tuples y^"^ leaves x invariant. This implies the following 
bound for Bn{X): 

n k—1 

B^{X)<nlY, H ^,5(ni > > l)e-^(»"-'^"=) . (E24) 

fc=i {yi,--- ,yk}cx {m,--- ,nfc)eP(fe,n) i=i 

For fixed I < k, each partition (ni, • • • ,nk) which contributes to the sum has to fulfill the condition 
I + k — 1 < n. The constrained given by the "delta-functions" 6{ni > 1) • • • (5(ri/_i > 1) can be replaced 
by the weaker constraint 5{l + k — 1 < n) which yields an upper bound of the right hand side: 

n fe-1 

B^{X) <nlJ2S{k,n)YS{l + k-l<n) ^ e-"^'^"' (£25) 



For each / < fc the sum over fc-elementary subsets can be estimated as follows: 



< J2 ^(^' k-l + 1, r)e-'^ 



(E26) 



where we have used Lemma 14 This yields the desired bound 

n k — 1 



Bn{X) <n\J2 Sik, n) ^ qik-i+i]Nk-i+i {X^/^ (E27) 



k=l 1=1 

which completes the proof of the lemma. 



Proof of Theorem^ By taking advantage of the expansion of Proposition |5] we can write the state ujx — 
Qjf^x £^gj. assume first, that the operators ai, • • • . a„ belong to the kernel of the single site 

restriction w. We first bound the correlation function of the functional Fx- 

|Fx(ai(8)---(g)a„)| 

|Ran(a:)|-l k~l 
a;eX" k=l 1=1 

|Ran(2;)|-l k-1 
kGX" k=l 1=1 

(E28) 

where we have used the fact that w(a^) = whenever the pre-image = 1 contains only one 



element. By inserting the definition of the quantities Bn{X) we find from Lemma 15 



<{Fx) < |X|-"/2GoB„(X) . (E29) 

By Proposition |6] the semi-norms i^^{uj'^^) are uniformly bounded, i.e. for each n e N, the constant 
An supxcA'^ni'^'^^) < oo is finite. Therefore we obtain for the induced state ujx the semi-norm 
bound: 

< A„ + GoB„ (E30) 
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Hence the semi-norm zy„ can be bounded by i/^ according to Lemma 1 1 



(E31) 



By applying this bound to ljx we obtain 



fc=0 



(E32) 



This implies that supj^^^ '^ni'^x) < oo and the net {cox)xcA has weakly ^n-fluctuations. 



